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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC),  under  Program  Element  65807F. 
The  results  were  obtained  by  ARO,  Inc.  (a  subsidiary  of  The  Sverdrup  Corporation), 
contract  operator  of  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee.  The  work  was 
done  under  ARO  Projects  No.  R33P-G0A  and  R33A-02A.  The  authors  of  this  report 
were  R.  C.  Bauer  and  J.  H.  Fox,  ARO,  Inc.  The  manuscript  (ARO  Control  No. 
ARO-ETF-TR-76-121)  was  submitted  for  publication  on  October  13,  1976. 
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1.0  INTRODUCTION 

The  afterbody  drag  of  a jet-propelled  vehicle  is  determined  by  the  strong  interaction 
of  the  engine  exhaust  with  the  external  stream.  Analytical  methods  to  compute  afterbody 
flow  fields  vary  considerably  depending  on  the  free-stream  Mach  number  and  whether 
the  external  flow  is  attached  or  separated.  For  conditions  with  attached  external  flow, 
Fig.  la,  an  "effective"  body  modeling  can  be  useful  although  there  is  an  unsolved  problem 
of  how  to  represent  the  jet  plume.  There  are  two  types  of  external  flow  separation:  (1) 
that  which  occurs  when  the  boundary  layer  cannot  sustain  the  adverse  pressure  gradient 
along  the  afterbody  surface,  Fig.  1 b,  and  (2)  that  caused  by  an  abrupt  change  in  afterbody 
geometry  (Fig.  lc).  The  second  type  of  flow  separation  is  more  amenable  to  analysis 
because  the  separation  point  is  fixed.  The  latter  separated  flow  is  of  practical  interest 
for  many  missile  configurations  which  operate  at  supersonic  Mach  numbers.  The  first  type 
of  flow  separation  is  difficult  to  analyze  because  of  inadequate  models  for  predicting 
boundary-layer  separation  phenomena,  particularly  on  aircraft  flying  at  transonic  and 
subsonic  Mach  numbers.  In  this  report,  only  flow  separation  caused  by  an  abrupt  change 
in  afterbody  geometry  will  be  considered  for  supersonic  external  stream  Mach  numbers. 
The  analysis  may  then  be  extended  for  application  to  the  first  type  of  separation  and 
for  all  Mach  numbers. 

The  basic  analytical  approach  followed  in  this  report  is  that  developed  by  Chapman 
and  Korst  (Ref.  1).  The  work  of  Korst  is  more  closely  followed  since  it  deals  with  turbulent 
base  flow.  Other  investigators  such  as  Addy,  Ref.  2,  and  Fong,  Ref.  3,  have  applied  the 
Chapman-Korst  theory  to  the  identical  problem  considered  herein.  All  subsequent 
investigators  have  recognized  that  the  major  discrepancy  in  the  Chapman-Korst  theory 
is  the  modeling  of  the  recompression  process.  Addy,  Ref.  2,  clearly  describes  how  the 
recompression  process  has  been  represented  by  the  various  investigators.  His  description 
is  quoted  in  the  following  two  paragraphs. 

"Since  the  recompression  process  is  not  well  understood  and  does  strongly  influence 
the  solution  to  the  base-flow  problem,  the  need  to  improve  agreement  between  predicted 
and  experimental  values  has  resulted  in  several  empirical  modifications  to  this  part  of 
the  flow  model.  Nash,  Ref.  4,  proposed  and  determined  an  empirical  "reattachment 
condition"  for  the  supersonic  single-stream  two-dimensional  case  which  was  approximately 
independent  of  the  supersonic  approach  Mach  number.  An  approach  was  adopted  by 
Carriere  and  Sirieix,  Ref.  5,  in  which  they  proposed  and  meticulously  determined  an 
empirical  law  of  reattachment,  for  a negligible  initial  boundary  layer,  based  on  an  "angular 
criterion  of  turbulent  reattachment."  Their  proposed  law  of  reattachment  is  shown  to 
correlate  well  the  data  for  reattachment  on  a wall  of  supersonic  axisymmetric, 
two-dimensional,  and  conical  flows.  Page,  Kessler,  and  Hill,  Ref.  6,  offer  an  alternative 
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empirical  reattachment  for  two-dimensional  supersonic  flow  based  on  the 
"discriminating-streamline"  velocity  ratio." 

"For  the  axisymmetric  two  stream  base- flow  problem  under  consideration  herein, 
the  recompression  resulting  from  the  interaction  of  separated  supersonic  free-stream  nozzle 
flows  apparently  cannot  be  analyzed  by  a direct  carry-over  of  the  empirical  reattachment 
criteria  determined  in  the  aforementioned  investigations.  As  a result,  Sirieix,  Delery,  and 
Mirande,  Ref.  7,  have  presented  a more  generalized  "criterion  of  turbulent  reattachment" 
which  involves  geometric  as  well  as  flow  variables.  This  criterion  is  the  basis  for  a limited 
comparison  between  predicted  and  experimental  values  of  the  axisymmetric  jet-on  base 
pressure  ratio.  Agreement  was  improved  for  the  cases  presented.  More  recently,  Dixon, 
Richardson,  and  Page,  Ref.  8,  considered  this  problem  under  the  experimentally  supported 
assumption  that  the  external  inviscid  flow  boundary  remains  straight  after  initially  turning 
to  adjust  to  the  base  pressure.  This  assumption  necessarily  results  in  an  axial-pressure 
gradient  from  which  the  internal  flow  boundary  is  determined.  The  corresponding  inviscid 
flow  fields  thus  determined  are  linked  to  an  approximate  axisymmetric  mixing  component 
by  a recompression  criterion  based  on  Goethert's  modification  of  Korst's  escape  criterion. 
This  approach  also  improves  the  agreement  between  the  predicted  and  experimental  results 
for  the  cases  presented  in  their  paper." 


a.  Attached  flow 


b.  Boundary-layer  separation  on  afterbody 
Figure  1.  Types  of  nozzle-afterbody  flows. 
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c.  Separation  because  of  abrupt  change  in  afterbody  geometry 
Figure  1.  Concluded. 


Addy,  Ref.  2,  presents  a relatively  simple  empirical  recompression  factor  for  cylindrical 
afterbodies  that  is  only  a function  of  the  ratio  of  afterbody  radius  to  jet  radius. 
Unfortunately.  Addy's  empirical  recompression  factor  had  to  be  modified  for  conical 
afterbodies  (Ref.  9).  Apparently,  a universal  empirical  recompression  factor  has  not  yet 
been  developed. 

In  this  report  a relatively  simple  recompression  flow  model  is  developed  which  yields 
a theoretical  value  of  the  recompression  factor.  The  various  factors  which  affect  the 
recompression  process  are  theoretically  approximated.  The  initial  boundary  layers  are 
included  in  the  turbulent  mixing  analysis  which  determines  the  development  with  distance 
of  the  mixing  process.  The  rate  of  mixing  is  based  on  turbulent  kinetic  energy  concepts 
and  accounts  for  all  differences  in  thermodynamic  properties  between  the  base  and  the 
corresponding  external  stream.  Finally,  the  analysis  includes  the  effects  of  a third  base 
bleed  gas. 
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2.0  THEORETICAL  DEVELOPMENT 

The  type  of  afterbody  flow  analyzed  is  that  in  which  the  flow  separates  because 
of  an  abrupt  change  in  afterbody  geometry.  A typical  streamline  pattern  associated  with 
this  type  of  afterbody  flow  is  shown  in  Fig.  lc.  In  Fig.  lc  the  streamlines  that  separate 
from  the  afterbody,  computed  by  assuming  no  mixing,  are  indicated  by  dotted  lines  and 
are  referred  to  hereinafter  as  the  inviscid  jet  boundaries.  The  turbulent  mixing  process 
shifts  the  inviscid  streamlines  to  the  solid  lines,  referred  to  as  the  dividing  streamlines, 
which  separate  the  high-speed  flow  from  the  entrained  flow.  The  intersection  of  the  two 
flow  fields  is  referred  to  as  the  recompression  zone  or  region.  In  the  recompression  zone 
a slip  line  is  formed  and  an  adverse  static  pressure  gradient  exists.  A portion  of  each 
turbulent  mixing  zone  will  not  have  sufficient  mechanical  energy  (total  pressure)  to 
penetrate  the  adverse  static  pressure  gradient  and  will  be  reversed  into  the  separated  region. 
The  stagnating  streamlines  are  those  streamlines  which  separate  the  flow  that  is  reversed  into 
the  separated  region  from  the  flow  that  has  sufficient  mechanical  energy  to  penetrate  the 
adverse  static  pressure  gradient  and  proceed  downstream.  For  typical  afterbody  flows,  as 
shown  in  Fig.  lc,  the  dividing  streamline  in  the  turbulent  mixing  zone  of  the  external  stream 
does  not  have  sufficient  mechanical  energy  to  proceed  downstream  and  is  reversed  into  the 
separated  region,  thus  forming  a vortex.  As  a result,  a portion  of  the  external  flow  is 
pumped  into  the  separated  region  and  must  be  pumped  out  by  the  turbulent  mixing  along 
the  jet  boundary.  To  do  this,  the  dividing  streamline  in  the  turbulent  mixing  zone  of  the  jet 
flow  must  have  sufficient  mechanical  energy  to  proceed  downstream.  Thus  the  equilibrium 
conditions  in  the  separated  region  are  those  which  make  the  mass  flow  pumped  in  equal  to 
the  mass  flow  pumped  out  of  the  separated  region.  In  addition,  conservation  of  energy  must 
simultaneously  be  established. 

The  application  of  the  Qiapman-Korst  base  pressure  theory  to  the  afterbody  flow 
shown  in  Fig.  lc  requires  a simplification  of  the  flow  field  shown  in  Fig.  2.  The  major 
simplifying  assumption  is  that  the  double  vortex  pattern  can  be  replaced  by  a region  of 
constant  fluid  properties  in  which  the  velocities  are  negligible.  In  addition,  the  base  region 
is  highly  turbulent  and  therefore  can  be  assumed  to  be  highly  stirred  and  thus  homogeneous. 
Based  on  this  assumption,  the  turbulent  mixing  processes  can  be  estimated  by  a 
single-stream,  constant  pressure  mixing  analysis.  Also,  the  constant  pressure  boundary 
condition  can  be  used  to  estimate  the  inviscid  portions  of  the  external  and  jet  flows. 
The  complete  flow  field  up  to  the  beginning  of  the  recompression  zone  can  then  be 
obtained  by  superimposing  the  mixing  zones  on  the  inviscid  flow  boundaries. 
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In  this  application  of  the  Chapman- Korst  theory  the  following  assumptions  are  used: 

1.  The  inviscid  portions  of  both  the  external  and  the  jet  flow  are  supersonic. 

2.  The  base  bleed  has  zero  momentum. 

3.  All  the  gases  obey  the  Perfect  Gas  Law. 

4.  All  the  gas  properties  are  constant  throughout  the  base  region  which  is 
defined  to  be  bounded  by  the  afterbody  wall,  the  lower  edges  of  the  mixing 
zones  (0  = 0.016),  and  the  beginning  of  the  recompression  region.  In 
addition,  the  velocities  in  this  region  are  negligible. 
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The  basic  equations  to  be  solved  are  those  for  the  conservation  of  mass  and  energy 
for  the  base  region  control  volume  shown  in  Fig.  2. 

The  equation  for  conservation  of  mass  is 

Yd  Yg  Y Y 

mE  - ^oo  J pudYj  + 2171^  J 1 pudYj  _ 2m00  J 2 pudY„  + 2 f 2 pudY  = 0 

(1) 


where  the  positive  sign  indicates  mass  flow  into  the  base  region.  This  equation  was  written 
assuming  the  beginning  of  the  recompression  is  at  the  point  (x,,.,  rj  in  order  to  simplify 
its  evaluation.  However,  this  does  not  introduce  any  error  since  Eq.  (1)  can  be  written 
as 

Yg  Yg 

mE  ~ 2nToc  X 1 Put^  i + 2nT«>  X Pu<^  2 = 0 (2) 

YDi  yl2 

Each  integral  term  in  Eq.  (2)  is  the  mass  flow  between  the  stagnating  and  dividing 
streamlines  which  must  be  the  same  for  the  mixing  zones  located  at  the  actual  beginning 
of  the  recompression  zone  (Fig.  2). 


The  equation  for  the  conservation  of  energy  in  the  base  region  is 

fVSi  Yg 

CpFTiEmE  + 2nr°°  J CpT>P“^  1 + 2ffr~  X 2 CpTtpudY2 


— 2 77  r „C_  T, 


yd  y 

X 1 pudY  + X Dz  P“dY2 

YT  y. 


= 0 


(3) 


Again,  the  mixing  zones  properties  are  evaluated  at  the  point  (x^,  rj  rather  than  at  the 
beginning  of  the  recompression.  This  assumption  introduces  compensating  errors  because 
the  first  two  integral  terms  overestimate  the  energy  into  the  base  region  by  the  same 
amount  the  last  two  integrals  overestimate  the  energy  leaving  the  base  region. 


The  evaluation  of  Eqs.  (2)  and  (3)  requires  complete  definition  of  the  inviscid-viscous 
flow  field  and  the  locations  of  the  dividing  and  stagnating  streamlines.  In  this  analysis 
the  inviscid  flow  fields,  the  turbulent  mixing  processes,  and  the  recompression  processes 
are  analytically  estimated  separately  and  combined  by  superposition  to  obtain  the 
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information  necessary  to  evaluate  Eqs.  (2)  and  (3).  The  base  region  properties  are  iterated 
until  Eqs.  (2)  and  (3)  are  satisfied. 

2.1  INVISCID  FLOW  FIELDS 

For  an  assumed  base  pressure,  Po  the  inviscid  flow  field  of  both  the  jet  and  external 
stream  are  computed  by  the  well-known  method  of  characteristics  solution  of  the  general 
potential  flow  equation  for  supersonic  flow.  Both  flows  are  assumed  to  separate  from 
the  physical  comers  of  the  afterbody  thus  neglecting  the  displacement  effect  of  the  initial 
boundary  layers.  The  major  information  to  be  obtained  from  this  calculation  is  the  location 
of  the  impingement  point  (x^  rj  of  the  inviscid  flow  boundaries,  the  length  of  each 
boundary  to  that  point,  and  the  inviscid  flow  conditions  along  the  high-speed  edge  of 
each  mixing  zone  after  superposition.  The  impingement  point  (x^,  rj  is  a point  on  the 
slip  line  that  is  formed  between  the  two  flow  fields.  The  slip-line  angle  and  maximum 
static  pressure  are  determined  by  the  interaction  of  the  inviscid  flows  along  the  high-speed 
edge  of  the  mixing  zones  at  (xm,  rj.  This  is  an  extension  of  the  method  developed  in 
Ref.  9 to  account  for  the  nonuniformity  of  the  inviscid  flow  on  the  recompression  process. 

2.2  TURBULENT  MIXING  ANALYSIS 

The  analysis  of  each  of  the  two  mixing  zones  is  based  on  the  following  assumptions: 

1.  The  mixing  process  is  isobaric. 

2.  The  width  of  the  mixing  zone  is  small  compared  to  the  radius  (b/r  < 0.3) 
so  that  the  locations,  within  the  mixing  layer,  of  the  dividing  streamline 
and  inviscid  boundary  are  the  same  as  for  two-dimensional  mixing  (Ref. 

10). 

3.  For  given  base  and  inviscid  flow  conditions,  the  eddy  viscosity  is  equal 
to  that  for  fully  developed  mixing. 

4.  The  inviscid  flow  is  uniform  in  the  vicinity  of  the  mixing  zone  and  has 
a velocity  equal  to  its  boundary  velocity. 

5.  The  Prandtl,  Lewis  and  Schmit  numbers  are  unity. 

2.2.1  Velocity  Profile 

The  velocity  profile  at  any  station  in  the  mixing  region  is  given  by  the  following 
equation  derived  by  Korst  in  Ref.  1. 
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4>  - - [i 
2 


erf  (77  - >jp)]  + — J 

V77 


1.0 


(4) 


where 


1^ 

n (initial  boundary  layer)  (5) 

Equation  (4)  defines  a parametric  family  of  velocity  profiles  in  terms  of  t/p  ranging  from 
the  initial  boundary-layer  profile,  for  t/p  = to  the  fully  developed  mixing  profile 
0 = 1/2(1  + erf  tj)  for  t/p  = 0.  An  analytical  method  has  been  developed  in  Ref.  11 
for  estimating  the  variation  of  rjp  with  (x/S).  However,  this  method  involves  a numerical 
integration  along  the  length  of  the  mixing  zone  and,  as  a result,  requires  significant 
computer  time.  A more  approximate  relationship  can  be  based  on  assumption  (3). 
For  two-dimensional  constant  pressure  mixing  the  relationship  is 


’Ip 


(6) 


where  a is  the  similarity  parameter  for  fully  developed  mixing.  A comparison  of  Eq. 
(6)  with  Chapman's  incompressible  experimental  data  (Ref.  12)  is  presented  in  Fig.  3a. 


a.  7/p  vs  x/8 

Figure  3.  Comparison  of  theoretical  tjp  and  0D  with  experiment. 
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The  theoretical  value  of  i?P  is  shown  always  to  be  greater  than  experiment  which  results 
in  an  under  prediction  of  the  velocity  ratio  on  the  dividing  streamline  as  shown  in  Fig. 
3b.  Included  in  Fig.  3b  is  the  "universal"  prediction  by  Alber  and  Lees  (Ref.  13)  which 
overpredicts  the  dividing  streamline  velocity  ratio  by  about  the  same  amount  the  present 
analysis  underpredicts  the  velocity  ratio.  Equations  (4)  and  (6)  are  considered  to  provide 
a reasonable  representation  of  the  velocity  profile  development  of  a turbulent  mixing  zone 
strongly  influenced  by  an  initial  boundary  layer. 


b.  <po  vs  x/5 
Figure  3.  Concluded. 

2.2.2  Thermodynamic  Composition 


The  thermodynamic  properties  of  interest  in  the  mixing  zone  are  total  temperature, 
Tt;  the  gas  constant  Rg,  and  the  sepcific  heat  for  constant  pressure  Cp.  As  a consequence 
of  assumption  5 the  thermodynamic  properties  at  any  point  in  the  mixing  zone  are 
determined  by  the  simple  mixture  equations  for  two  gases,  i.e.,  the  high-speed  gas  and 
the  base  region  gas.  The  equations  are 


and 


and 


k + (1 


k) 


R 


f=oo 


k + (1 


k) 


(7) 


(8) 


r, 


k x (i  _ it) 


(9) 
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where  the  subscript  «.  indicates  the  high-speed  stream  and  k is  the  local  mass  fraction 
of  the  high-speed  stream  at  any  point  in  the  mixing  zone. 


The  distribution  of  k through  the  mixing  zone  can  be  determined  by  simultaneously 
solving  the  conservation  of  mass  and  momentum  equations.  The  conservation  equations 
for  a typical  mixing  zone,  shown  in  Fig.  4,  are  Conservation  of  Mass  Flow: 


Figure  4.  General  mixing  zone. 


(10) 


(11) 


(12) 
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Since 


and  for  Y < Y^ 


y = Y _ Y 

* 30 

dy  = dY 

Piui  = Pi“i2  = o 


for  Y > Ym  + 5 


Pi“i  = P» 


PiUi  = PooUoo 


Equation  (12)  can  be  written  as 


Therefore, 

k = PiUi  _ PiUi2 

pu  pUUoo 

The  bounds  on  PjUi  also  apply  to  Eq.  (14),  i.e., 


for  Y < Y.. 


for  Y_  < Y < (Y.  + 5) 

k = PiUi  _ PiUi2  u 

Pu  Puu<»  u=o 


(13) 


(14) 


(15) 


(16) 


and  for  Y > (Y^  + 6) 


(17) 


It  should  be  noted  that  for  zero  initial  boundary  layer  k = u/u^,  which  is  the  well-known 
Crocco  relation  for  unity  Prandtl,  Lewis  and  Schmit  numbers. 

In  terms  of  the  nondimensional  mixing  variables. 


15 


AEDC-TR-76-158 


Eq.  (11)  becomes 

*•  - ,u  - ,p  [t  - f]  - (?  - c-!) 

and  Eqs.  (15),  (16),  and  (17)  become 
For  ??  < i?oo 


k = d> 


(19) 


for  Voo  < V < <7oo  + V?) 


k = <f>  ~ 


0 

-Tt  /C 

K v 

] 

(1?- 

- 1 1 - 

n . 7n 

VP  - (V-riJ 

" 2 2" 
vpa  - cM2  (17-5700)  “ 

(20) 


for  i?  > Ci7oo  + Vp) 


k = 4> 


(21) 


Equations  (18)  through  (21)  are  used  to  determine  t}.,  the  location  of  the  inviscid  flow 
boundary  in  the  turbulent  mixing  zone,  which  is  used  to  properly  superimpose  the  mixing 
zone  on  the  inviscid  flow  field  and  to  establish  the  mass  fraction,  k,  distribution  through 
the  mixing  zone. 


2.2.3  Location  of  Dividing  Streamline 

The  location  of  the  dividing  streamline  in  the  mixing  zone  is  obtained  by  applying 
the  conservation  of  mass  condition  to  the  mixing  process  shown  in  Fig.  4.  The  basic 
equation  is 

8 v,j 

J PjUidy  + Poouco  tYU  - Y«-  = J pud)'  (22) 

0 yd 

Since  Eq.  (22)  involves  it  is  necessary  to  combine  Eq.  (22)  with  the  conservation 
of  momentum  equation  [Eq.  (1 1)]  to  yield 

8 8 Yy  Yy 

uoo  f Pi«i«ly  - J Piui2 dy  + I pu2dY  = !!„  J pudY  (23) 

0 o yl  yd 
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By  definition 


[S  - 5*]  = J P;U; 


dy 


and 


Pmu»o2[S  - S*  - #1  = J PiUj2dy 

o 

Thus,  Eq.  (23)  becomes 

Yu  Yu  yd 

P<*uoo  6 + J pu^d\  = J"  pudY  - ijro  j pudY 


In  nondimensional  form  Eq.  (26)  is 

"i’id  '"-'iu-  - TT7EJ,  (1) 

and  Eqs.  (24)  and  (25)  become 

T=1-(1-  C~H  TTEJtJ  A(j) 

and 

l.d  Jp) 

« o 1 - C^^.2  \5J 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


2.2.4  The  Similarity  Parameter,  a 

In  the  previous  section  the  mixing  velocity  profile  is  defined  in  terms  of  the  position 
parameter,  ijp.  However,  to  determine  7jp,  the  similarity  parameter,  o,  for  planar,  fully 
developed,  constant  pressure  turbulent  mixing  is  determined  by  applying  the  turbulent 
kinetic  energy  (TKE)  method  to  the  mixing  zone  shown  in  Fig.  4 for  no  initial  boundary 
layer.  The  planar  value  is  then  corrected  for  axisymmetric  effects.  The  basic  equation 
is 

i A A A 

—j—  <YU  - Y«>  - T K pu  dY  = 3V  puKdY  4 J XDdY  (30) 

z H ' L yi. 
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where 


(u)2 


(u02  + (u)2  + 


(31) 


and  D = dissipation  per  unit  volume 


Equation  (30)  states  that  the  loss  in  mean  flow  kinetic  energy  (production  of  TKE)  equals 
the  convection  of  TKE  at  station  x plus  the  total  dissipation  in  the  region  between  x 
= 0 and  x.  For  isotropic  turbulence  the  dissipation  is 


D = 


(32) 


By  definition 


(33) 

(34) 


Substituting  into  Eq.  (30)  yields 


Po cuo 


1 J7U  Jv  /a 0a\ 

- V«)  — 2 5 p"3d,?  = I puKd,?  + / f -=— ) PK3/2dl?  (35) 

»?L  *7l  11  \Vh/ 


Dividing  by  pjO 


■ t C (f>  • C Q'O*  • C®  (£)&'  * 


(36) 


For  a perfect  gas 


P = — . = 


#«T  R*T'(1 


rpj 


and 


P«> 


Boc 


Pc 

T,  Cl  - C J\ 

oo 
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U - c,2) 

(1  - Ca) 


By  definition 


and 


2CpTt 


C 


2 

OO 


2 Cr 


C 


2 


or 


(37) 


(38) 


(39) 


Substituting  Eq.  (39)  into  Eq.  (36)  yields 


Substituting  Eq.  (18)  for  no  boundary  layer  into  Eq.  (40)  yields 
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Equation  (41)  is  the  TKE  equation  which  is  solved  by  using  the  empirical  information 
for  a2  and  K developed  by  Peters  (Ref.  14). 


Let 

Y K 

= XT  (42> 

where 

Substituting  Eq.  (42)  into  Eq.  (41)  yields 


Assume: 


or 


rm  = alPmKm  (sflear  stress  at  4>  = -j-) 
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also  assume  ai,  a2,  a and  are  independent  of  17. 
Therefore,  Eq.  (43)  becomes  as  follows 


(45) 


The  shear  stress  distribution  in  the  mixing  zone  can  be  determined  by  applying  the 
Conservation  of  Momentum  Equation  as  in  Ref.  1. 


The  general  equation  is 


- * to,),  - 


Poo^o 


Therefore,  at  <f>  - 1/2  and  17  = 0 for  fully  developed  mixing, 


<1-0  1 
— [<I2)0  - - [<I,>0  - (1,)^]] 


Poc*a 


Substituting  Eq.  (46)  into  Eq.  (45)  yields 


(46) 


OH  -• 

For  the  present  set  of  variables,  Peters'  distribution  of  K in  the  mixing  zone  is 
For  r,  < - K = o 


(47) 
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K = 0 for  n > 

2 

The  two  remaining  unknowns  in  Eq.  (47)  are  ai  and  a2 . Peters  assumes  ai  has  a universal 
value  of  0.3,  and  a2  is  a function  of  the  turbulent  Reynolds  number  defined  as 


u„  b , _ 

Rt  = (50) 

€m 

Since 

du 

t = p e — 
m r m m 

em  in  Eq.  (50)  can  be  eliminated  to  yield 


Substituting  Eqs.  (33)  and  (46)  into  Eq.  (52)  yields 


(51) 

m 
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For  fully  developed  turbulent  mixing  the  velocity  profile  assumed  by  Peters  is 


*-T 


lb, 


i-R^ost  - nf- 

2 L L2  V 

dd>  I tt  / I \ 

drJ  In.  Zrlh  \ 2 / 


From  Eqs,  (39),  (7),  (8),  and  (9) 


l + _ — (l-r  *) 


' V 

l + ■ — — 


R 


±L  W c-! 

2rTO“ 


Substituting  Eqs.  (55)  and  (56)  into  Eq.  (53)  yields 


(54) 


(55) 


(56) 


2 TT<J  ( 1 - • 


R-y  = 


-m 


C_  T, 


c rr, 

Poo  ^oaJ 


- c. 


(I2>0 


Peters  empirically  defines  a2  as 


For  o < rt  < 70 


For  70  < Rt  < 145.223 


For  R'j-  > 145.223 


0.99  + 0.0]  Rt 


1.69 


a.,  = 


3.3  — 5.  ]9£-0  °0806Rt 


12  C 

where  c is  empirically  defined  as 


7D 


(57) 


(58) 


(59) 


(60) 
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For 


U - C*,2)  > 1.0 


c 


0.984  + 0.016 


(61) 


For 


R*.T«. 


°c 


(1  - O < 1.0 


c 


0.95  + 0.05 


(62) 


The  limits  ijl  and  r?u  on  the  integrals  Ij,  I2 , I3 , and  I4  are  taken  to  be  -3.0  and  +3.0, 
respectively,  to  obtain  an  accuracy  equal  to  that  obtained  by  Korst  (Ref.  1).  These  limits 
correspond  to  velocity  ratios  of  0.00001105  and  0.99998895,  respectively,  for  the  fully 
developed  mixing  velocity  given  by  0 = 1/2(1  + erf  1?).  Peters  used  the  following  equation 
to  determine  o: 

3.12 

~db~  (63) 

dx 

From  Eq.  (33) 

db 

17  = T <64> 

Substituting  Eq.  (64)  into  Eq.  (63)  yields 

Vb  = 3-12  (65) 

Based  on  these  values  for  17L » iju>  and  rj b the  variation  of  a with  Mach  number  computed 
by  this  method  is  compared  with  that  by  Peters  in  Fig.  4.  As  shown,  a determined  by 
the  present  method  is  about  ten  percent  less  than  that  determined  by  Peters.  This  is  believed 
to  be  caused  by  the  use  of  an  error  function  velocity  profile  rather  than  the  cosine  profile 
used  by  Peters.  Included  in  Fig.  5 is  the  estimated  a variation  with  Mach  number  based  on 
the  well-known  Donaldson  and  Gray  correlation  of  eddy  viscosity,  Ref.  21.  The  TKE 
method  was  used  in  the  present  afterbody  analysis  because  the  empirical  constants  are 
believed  to  be  more  general  than  those  based  on  the  eddy  viscosity  concept. 
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Figure  5.  Variation  of  similarity  parameter  with  free-stream  Mach  No. 


2.2.5  Axisymmetric  Effect  on  a 

Although  the  assumption  that  the  width  of  the  mixing  zone  is  small  compared  to 
the  radius  does  not  eliminate  an  axisymmetric  effect  on  the  similarity  parameter,  a,  it 
does  allow  a simple  correction  to  be  applied.  The  correction  factor  is  derived  for  the 
mixing  along  the  jet  plume  boundary.  The  jet  plume  boundary,  indicated  in  Fig.  6,  is 
assumed  to  be  a conical  surface.  The  condition  applied  is  that  the  momentum  of  the 
entrained  mass  flow  equals  the  total  shear  force  along  the  dividing  streamline.  The  basic 
equation  is 

— cos  fi  = JY  pu2  [r  + [r  - (V  - YJ  cos  jQ]]  dY  (66) 

Let 


v 


O * ^ <:<>  s ft 

— 

X 
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where 


or 

Since  by 


axisymmetric  o 

rxCrjj  + r) 


/3  ff  A cos  £ r)L 


J 

2,-(- 

-)  (*?  - 

Vl 

L V 

V 

drj 


JV]}  / x \ f^D  2 

pu2ir)  - ( ) J pn  ~ dij 

*)L  W/  V 


definition 


cta1>  cos  /3 
»?b  = 5 


Eq.  (68)  can  be  written 
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By  assumption  2 
For  b/r  < 0.3 

(t)  ^ pu2(  ~ 7?o°)d’7  <<:  ^ pt,2d'7 

The  condition  expressed  by  Eq.  (70)  is  numerically  verified  by  the  following  example: 
For  — = 0-3 

r 

Moo  - 3.0 

From  Ref.  10 

rjjj  = 0.3115 

Voo  = 0 7388 

’Jd 

„ 2j_  0.23 

P u a>?  = — 

-3-0  Poo%c 

Vo 

| pu2rfdri  = 0 
J-3 . o 

and  also  assume 

/3  = 0° 

% = 3.12 

Substituting  into  relation  (70)  gives 

-0.0335  « 1.0 

thus  verifying  relation  (70)  for  this  set  of  conditions. 

Therefore,  Eq.  (69)  can  be  written 


raA 


2 


(71) 
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If  the  mixing  were  two-dimensional  the  basic  equation  would  be 


and  in  the  transformed  plane 


(72) 


(73) 


Substituting  Eq.  (73)  into  Eq.  (71) 


2 


TO 


Letting  R = -7- 


yields 


(74) 


In  a similar  manner  the  correction  factor  for  the  external  stream  mixing  can  be  shown 
to  be 


■ (ifh) ' <75> 

Equation  (74)  shows  that  a for  mixing  along  an  expanding  plume  is  always  greater  than 
that  for  two-dimensional  mixing,  whereas  for  a converging  plume  Eq.  (75)  shows  that 
a is  always  less  than  the  two-dimensional  value. 

2.3  RECOMPRESSION  ANALYSIS 

The  most  important  aspect  of  the  recompression  process  that  must  be  analytically 
modeled  is  that  the  magnitude  of  the  pressure  at  the  stagnation  point  is  between  the 
base  pressure,  pc,  and  the  maximum  recompression  pressure,  p3,  is  illustrated  in  Fig.  7 
for  the  case  of  rearward-facing  step.  Originally,  Chapman  and  Korst  assumed  the  stagnation 
pressure  to  be  the  peak  pressure,  p3 , but  this  was  disproved  by  Nash  (Ref.  4).  Recently 
two  observations  have  been  published  concerning  the  recompression  process  that  suggests 
the  phenomena  can  be  realistically  represented  by  a relatively  simple  analytical  model. 
The  two  observations  are: 
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Distance  along  Wall,  x 


Figure  7.  Typical  recompression  process. 


1.  Narayanan,  et  al..  Ref.  15,  experimentally  show  that  for  incompressible  planar 
flow  the  pressure  distribution  in  the  recompression  zone  has  "almost  perfect  similarity" 
using  the  similarity  variables  (p  - pc)/(P3  - Pc)  and  (x  - x*)/h.  In  the  present  analysis 
it  is  convenient  to  use  the  similarity  variable  (x  - Xi)/(X2  - xj)  which  is  related  to  (x 
- x*)/h  by 


(76) 
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From  the  experiments  presented  in  Ref.  15 


Fig.  8 compares  the  experimental  data  of  Ref.  15  with  the  assumed  relation  between 
the  similarity  variables. 


(77) 


x-Experiments  (Ref.  15) 
h = Step  Height 

P ' Pc 


x = x for 


= sin1 


P3’Pc 


1 

2 

+ 2 


Figure  8.  Illustration  of  similarity  in  the  recompression 
pressure  distribution. 


Since  values  of  the  similarity  variables  (p  - pc)/(p3  - pc)  and  (x  - x!)/(x2  - Xi)  are  always 
between  0 and  1.0  it  is  reasonable  to  assume  Eq.  (77)  also  applies  for  planar  supersonic 
flow.  This  is  verified  by  the  theoretical  and  experimental  pressure  distributions  presented 
by  Alber  and  Lees,  Ref.  13,  and  by  experimental  results  of  Viswanath  and  Narasimha, 
Ref.  16. 


2.  Chow's  analysis  of  the  recompression  process  for  flow  over  a backward-facing  step 
presented  in  Ref.  17  shows  that  the  total  pressure  on  the  stagnating  streamline  continues 
to  increase  at  about  the  same  rate  all  the  way  to  the  stagnation  point  on  the  wall.  This 
result  is  presented  in  Fig.  9.  The  same  phenomena  was  reported  by  O'Leary  and  Mueller 
in  Ref.  18  from  a numerical  solution  of  the  Navier-Stokes  equation  for  incompressible 
laminar  flow  over  a backward-facing  step.  Based  on  these  results  it  is  assumed  that  the 
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Figure  9.  Estimated  properties  along  stagnating  streamline 
by  Chow  (Ref.  17). 


adverse  pressure  gradient  along  the  stagnating  streamline  in  the  recompression  zone  does 
not  significantly  alter  the  gradient  of  total  pressure  compared  to  a constant  pressure  shear 
layer.  Therefore,  the  total-pressure  development  through  the  recompression  zone  to  the 
stagnation  point  can  be  estimated  by  the  constant  pressure  mixing  analysis' neglecting  the 
adverse  pressure  gradient. 

2.3.1  Analytical  Model 

The  recompression  flow  field  is  modeled  in  a manner  similar  to  that  presented  in 
Ref.  9.  The  assumed  flow  field  is  presented  in  Fig.  10.  The  mixing  is  superimposed  on 
the  inviscid  flow  boundary  in  its  proper  position  as  determined  by  conservation  of 
momentum,  Eq.  (18).  The  reverse  flow  is  assumed  to  turn  symmetrically  about  a line 
emanating  from  the  stagnation  point  at  half  the  angle  of  the  stagnating  streamline.  The 
intersection  of  the  symmetry  line  with  the  lower  edge  of  the  mixing  zone  (defined  to 
be  where  <f>  = 0.016)  determines  the  beginning  of  the  recompression  process,  xi.  The 
end  of  the  recompression  region,  X2,is  defined  to  be  the  location  of  the  maximum  slip 
line  static  pressure,  P3 . Since  an  isobaric  turbulent  mixing  process  simply  redistributes 
the  momentum  of  the  inviscid  flow,  the  total-pressure  force  between  xj  and  X2  must 
equal  that  for  inviscid  flow.  Applying  this  condition  for  uniform  inviscid  flow  and  the 
slip  line  pressure  distribution  given  by  Eq.  (77)  yields  the  following  equation  for  x2 
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*2  = 2x«  ~ X1 


(78) 


Figure  11.  Conditions  jit  stagnation  point. 


The  analysis  was  applied  to  the  nozzle-afterbody  flow  shown  in  Fig.  2 based  on 
the  assumption  that  the  slip  line,  Fig.  2,  is  a conical  surface.  The  following  equations 
can  be  derived  from  the  geometry  of  the  jet  plume. 
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^1  {rloc]  cos“i  /I  + Ro 


= X..  _ ■ ■ ■ -1  /1  + R~\ 

S!  °°  O j sin  [fil  -dj-w]  \ 2R,*,  / 


(79) 


X j = X ^ - 

*1  S1 


cos  a i cos 


/x-  - xsN 

:oscj  — ^ — ^cos  (j8]  —cj]) 
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(80) 


x2i  . «_  - X 


where 


tan 


/ 


A~i-7sA  /,iR~\  jf-h 


(81) 


(82) 


The'- pressure  at  the  stagnation  point  is  given  by  the  following  equation  from  Eq.  (77), 


Pc 

Equations  (79)  through  (83)  are  evaluated  using  the  properties  of  the  mixing  zone  located 
at  (Xj,  R,).  The  only  unknown  is  tjs  i > the  location  of  the  stagnation  point  in  terms  of 
the  mixing  coordinates.  For  an  assumed  value  of  77s  A Eqs.  (79)  through  (83)  determine 
the  pressure,  p§  l , at  the  corresponding  stagnation  point.  Assuming  p§  j is  the  total  pressure 
of  the  streamline  impinging  at  r)slt  the  velocity  ratio,  fej,  can  be  determined  from 


- 1 


Pc 


sin 


xs,-  X, 


X,  -X, 


(83) 
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The  values  of  Tjg  ^ /T(  ^ , Cp  g ^ / Cp  ^ , Rgj/Rgj  are  functions  of  <j)  and  t?sj  defined  by  Eqs. 
(7),  (8),  and  (9). 


The  velocity  ratio,  <j>s1,  can  also  be  determined  from  Eq.  (9)  in  the  form 


**i-  («:)”'  <«) 


The  stagnation  point  location  is  defined  when  the  assumed  value  of  17 s j results  in  equal 
values  of  0s  j from  Eqs.  (84)  and  (86). 

The  recompression  of  the  external  stream  along  the  conical  slip  line  (Fig.  2)  is 
computed  independent  of  the  recompression  of  the  jet  plume.  As  a result,  the  location 
of  stagnation  point  of  the  external  stream  may  be  different,  as  may  be  the  stagnation 
pressure,  from  that  of  jet  plume.  Physically,  for  this  to  be  true,  it  would  be  necessary 
for  a secondary  vortex  system  to  exist  in  the  recompression  zone  between  the  two 
stagnation  points.  Since  the  existence  or  nonexistence  of  such  a vortex  system  has  not 
been  established,  the  authors  of  this  report  have  elected  not  to  impose  the  single  stagnation 
point  condition. 


The  recompression  equations  for  the  external  stream  are 

e2(7«.9  ~ ,?50)  c0sC0  cos  a2  /R2+r» 


Xc  = - 
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where 


(91) 


(92) 


(93) 


where 


where 


7p 


2 


(97) 


(98) 
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2.4  SOLUTION  CONDITIONS 


The  mixing  and  recompression  analyses  are  evaluated  based  on  values  of  pc,  Tc,  Rgc, 
and  CPc.  However,  it  is  not  necessary  to  assume  individual  values  for  each  of  the 
thermodynamic  parameters  since  they  are  determined  by  the  mixture  equations  for  the 
three  gases  entering  the  base  region.  These  equations  are: 


T.  Eq.  (100) 
ll 

These  equations  utilize  the  identity 


(99) 


(100) 


(101) 


(102) 


to  eliminate  kC2. 


+ 


1.0 


(103) 


Using  Eqs.  (99)  through  (102)  requires  assuming  the  three  variables  pc,  kci  and  kCE 
rather  than  the  four  variables  pc,  Tc,  Rc,  and  CPc.  In  addition,  for  most  afterbody  flows 
there  is  no  bleed  flow  so  that  only  two  variables,  pc  and  kc  j need  to  be  specified.  This 
greatly  simplifies  the  iteration  procedure. 


For  the  assumed  variables  to  be  a solution  the  conservation  of  mass,  Eq.  (1),  must 
be  satisfied  and  the  assumed  values  of  kC[  and  kC£  must  be  equal  to  the  values  computed 
from  the  equations  for  conservation  of  species. 


The  conservation  of  mass  equation  in  the  transformed  plane  is 


(104) 
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where 
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The  equations  for  kci  and  kCE  are  derived  as  follows: 


From  Fig.  2: 


(110) 


and 


m E 


CE 


1 ^ I C 
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In  the  transformed  plane,  Eqs.  (110)  and  (111)  are 


(111) 
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(112) 
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(113) 


3.0  EVALUATION  OF  THEORY 

The  purpose  of  this  evaluation  is  to  establish  the  general  quantitative  accuracy  of 
the  analytical  technique.  Evaluation  of  the  various  components  of  the  analysis  is  not 
possible  at  this  time,  since  the  existing  experimental  data  consist  of  only  base  pressure 
measurements.  The  experimental  data  selected  for  this  evaluation  are  those  which  show 
the  effects  on  base  pressure  of  (1)  base  bleed,  (2)  thermodynamic  properties,  (3)  initial 
boundary  layers,  and  (4)  geometry. 

In  the  following  sections  the  analysis  is  compared  to  experiments  having  an  external 
stream  Mach  number  of  2.0.  The  jet  Mach  number  is  also  2.0  except  for  the  one  hot 
rocket  configuration.  These  data  were  selected  because  the  present  analysis  does  not  include 
a method  for  estimating  the  change  in  initial  boundary-layer  properties  caused  by  the 
expansion  process  at  the  separation  point.  Many  approximate  analyses  have  been  developed 
for  this  phenomena,  Ref.  4,  and  all  show  that  at  Mach  2.0  the  momentum  thickness 
of  the  boundary  layer  is  relatively  unaffected  by  the  expansion  process. 

3.1  DAVIS  EXPERIMENT 

Davis,  Ref.  19,  investigated  the  effect  of  chemical  reactions  on  base  pressure  for 
M = 2.0  flow  over  a two-dimensional  blunt  base  by  injecting  hydrogen  in  the  base  region 
and  igniting  with  a pilot  flame.  Of  interest  in  this  report  are  the  data  presented  in  Ref. 
19  showing  the  effect  on  base  pressure  of  a hydrogen  bleed  without  chemical  reactions. 
Two-dimensional  base  flow  can  be  simulated  by  an  axis ym metric  configuration  having  a 
small  base  height  fa  - r^),  Fig.  2,  relative  to  the  jet  radius,  rjj.  Experimental  data  for 
an  axisymmetric  system  having  r2  - rN/rN  = 1.25  obtained  by  Reid  and  Hastings,  Ref. 
20,  are  shown  in  Fig.  12a  to  agree  well  with  the  Davis  two-dimensional  value  for  no 
bleed.  The  axisymmetric  configuration  selected  for  simulating  two-dimensional  base  flow 
has  a value  of  r2  - in/Tn  = 1.2. 


38 


AE  DC-TR-76-1 58 


P(D'  Ma,iH 


a.  Base  pressure 

Figure  12.  Comparison  of  theory  with  Davis  experiment. 

Davis  presents  his  data  in  terms  of  a base  bleed  parameter  defined  as 


II 


m 
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t,  <y„-i> 

DO 

2yx 


(114) 


The  analysis  in  this  report  requires  as  input,  kg,  the  ratio  of  bleed  mass  flow  to  jet  mass 
flow.  The  two  bleed  flow  parameters  are  related  by 


h 


(115) 


For  the  specific  conditions  of  the  Davis  experiment 

kp;  = 2.823  H (116) 

The  theoretical  variation  of  base  pressure  with  bleed  flow  is  compared  with  experiment 
in  Fig.  12a.  The  theory  agrees  well  with  experiment  for  no  bleed,  but  is  greater  than 
experiment  by  about  50  percent  for  the  maximum  bleed  case.  This  indicates  the  theoretical 
rate  of  mixing  is  too  low,  as  was  previously  shown  for  incompressible  mixing  in  Figs. 
3a  and  3b. 
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The  theoretical  variation  of  specific  heat  ratio  in  the  base  region  is  presented  in  Fig. 
12b  and  shows  the  expected  trend  with  bleed  flow.  The  variation  of  the  similarity 
parameter,  a,  with  bleed  flow  is  shown  in  Fig.  12c  to  decrease  with  increasing  bleed  flow. 


b.  Base  thermodynamic  properties 


c.  Mixing  similarity  parameter 
Figure  12.  Continued. 


The  trend  of  a decreasing  with  H2  bleed  flow  is  inconsistent  with  both  the  Donaldson 
and  Gray  and  Channapragada  formulations  for  turbulent  mixing  (Refs.  21  and  22).  Finally, 
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the  variation  of  a Nash-type  recompression  factor,  N,  Ref.  4,  with  bleed  flow  is  presented 
in  Fig.  12d.  For  no  bleed  flow  the  theoretical  value  of  0.42  agrees  well  with  the  value 
0.35  assumed  by  Nash  based  on  experiment.  The  decrease  in  the  recompression  factor, 
N,  with  increasing  bleed  flow  is  consistent  with  the  experimental  results  obtained  by 
Carriere  and  Sirieix,  Ref.  23,  using  air  as  the  bleed  gas. 


d.  Recompression  factor 
Figure  12.  Concluded. 


3.2  TWO-DIMENSIONAL  BACKWARD-FACING  STEP  EXPERIMENTS 

The  theoretical  effect  of  initial  boundary-layer  momentum  thickness  on  the  base 
pressure  for  Mach  2.0  flow  over  a backward-facing  step  is  compared  with  experiment  in 
Fig.  13a.  The  theoretical  results  were  obtained  for  the  same  axisymmetric  configuration 
used  to  represent  the  Davis  two-dimensional  experiment.  The  theoretical  base  pressure 
is  shown  in  Fig.  13a  to  agree  well  with  Chow's  analysis  (Ref.  17)  for  relatively  thin  initial 
boundary  layers,  and  with  Alber  and  Lees  analysis  for  the  thicker  initial  boundary  layers. 
The  theoritical  base  pressure  is  about  25  percent  larger  than  experiment  at  the  higher 
values  of  6J h.  The  theoretical  mixing  similarity  parameter,  o,  is  presented  in  Fig.  13b 
as  a function  of  boundary-layer  thickness,  but  its  variation  is  due  entirely  to  the  change 
with  base  pressure  of  the  inviscid  flow  Mach  number  along  the  edge  of  the  mixing  zone. 
The  theoretical  variation  of  a Nash-type  recompression  factor  with  initial  boundary-layer 
thickness  is  presented  in  Fig.  13c.  The  magnitude  of  the  recompression  factor  is  in  good 
agreement  with  that  presented  by  Nash  in  Ref.  4.  In  addition,  the  small  variation  of  the 
recompression  factor  with  boundary-layer  thickness  is  consistent  with  the  assumption,  by 
Nash,  of  a constant  recompression  factor  used  in  his  analysis. 
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a.  Initial  boundary-layar  effect  on  base  pressure 


b.  Mixing  similarity  parameter 
Figure  13.  Comparison  of  theory  with  two-dimensional 
backstop  experiment. 
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c.  Recompression  factor 
Figure  13.  Concluded. 


3.3  ROCKET-AFTERBODY  CONFIGURATION 

The  theory  was  applied  to  the  rocket-afterbody  configuration  shown  in  Fig.  14a  to 
determine  the  effect  on  base  properties  of  the  external  stream  boundary-layer  thickness. 


M2  = 2.05 


a.  Base  pressure 

Figure  14.  Application  of  theory  to  a rocket-afterbody  configuration. 
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The  theoretical  variation  in  the  base  pressure  shown  in  Fig.  14a  is  much  less  than  that 
computed  for  the  flow  over  a backward-facing  step,  Fig.  13a.  Apparently  the  base  heating 
effect  for  a thin  boundary  layer,  shown  in  Fig.  14b,  prevents  the  base  pressure  from 
decreasing  at  low  values  of  82. 

The  variation  of  the  mixing  similarity  parameter,  a,  and  the  recompression  factors 
for  each  stream  are  presented  in  Figs.  14c  and  d.  The  recompression  factors  are  different 


b.  Base  thermodynamic  properties 


c.  Mixing  similarity  parameters 
Figure  14.  Continued. 
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External  Flow  Boundary-Layer  Thickness  Ratio, 

d.  Recompression  factors 
Figure  14.  Concluded. 

because  they  are  computed  independent  of  each  other  as  previously  discussed.  The 
recompression  factor,  N,  indicates  the  relative  location  of  the  stagnating  streamline  in  the 
mixing  zone.  For  N = 1.0  the  stagnating  streamline  is  located  in  the  high-speed  portion 
of  the  mixing  zone.  For  N = 0 the  stagnating  streamline  is  located  at  the  low-speed  edge 
of  the  mixing  zone  and  the  entire  mixing  zone  mass  flow  will  proceed  downstream.  The 
recompression  factor  for  the  external  stream,  N2,  is  shown  in  Fig.  14d  to  increase  as 
the  boundary-layer  thickness  increases,  thus  reflecting  the  rejection  of  more  and  more 
of  the  boundary-layer  mass  flow  into  the  base  region.  As  a result,  the  recompression  factor 
for  the  jet,  Ni,  is  decreasing  indicating  an  increase  in  the  jet  pumping  capacity. 

3.4  REID  AND  HASTINGS  EXPERIMENT 

Reid  and  Hastings,  Ref.  20,  experimentally  studies  a family  of  afterbody  geometries, 
all  at  an  external  stream  and  jet  Mach  number  of  2.0.  The  most  extreme  configuration 
tested  had  an  external  surface  radius  five  times  the  jet  radius.  The  theoretical  base  pressure 
as  a function  of  external  boundary-layer  thickness  for  that  configuration  is  presented  in 
Fig.  1 5a.  The  theoretical  base  pressure  increases  significantly  with  boundary-layer  thickness 
just  as  it  does  for  flow  over  a backward-facing  step.  The  theory  is  about  35  percent  greater 
than  the  experimental  value  given  by  Reid  and  Hastings.  The  variation  of  the  recompression 
factor  for  each  stream  is  presented  in  Fig.  15b,  and  the  trend  is  the  same  as  for  the 
rocket-afterbody  configuration  previously  discussed. 
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External  Flow  Boundary-Layer  Thickness  Ratio, 


a.  Base  pressure 


External  Flow  Boundary-Layer  Thickness  Ratio,  62^1 
b.  Recompression  factors 

Figure  15.  Comparison  of  theory  with  Reid  and  Hastings  experiment. 

4.0  CONCLUDING  REMARKS 

The  evaluation  of  the  analysis  shows  that  usually  the  theoretical  base  pressure  is 
greater  than  experiment.  The  experimental  data  are  too  limited  to  indicate  which  parts 
of  the  analysis  are  most  seriously  in  error.  Since  the  only  empiricism  used  is  in  the  turbulent 
mixing  analysis,  this  would  be  the  simplest  to  modify.  However,  a higher  rate  of  mixing 
is  required  to  reduce  the  theoretical  base  pressure  which  could  be  obtained  by  use  of 
the  Donaldson  and  Gray  formulation  (Ref.  21)  for  the  turbulent  mixing  analysis. 
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APPENDIX  A 

DEFINITION  OF  INTEGRALS 
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NOMENCLATURE 

A Area 

A*  Throat  area  of  jet  nozzle 

ai  Turbulent  kinetic  energy  constant,  Eq.  (44) 

a2  Turbulent  kinetic  energy  constant,  Eq.  (34) 

b Width  of  mixing  zone 

C Crocco  number,  Eq.  (38) 

Cp  Specific  heat  at  constant  pressure 

c Turbulent  kinetic  energy  constant,  Eqs.  (61)  and  (62) 

D Turbulent  kinetic  energy  dissipation  per  unit  volume,  Eq.  (32) 

G Mass  flow  between  dividing  and  stagnation  streamlines 

H Base  bleed  flow  parameter,  Eq.  (114) 

h Step  height 

1 Integrals  defined  in  Appendix  A 

AI  Nondimensional  mass  flow  between  dividing  and  stagnating  streamlines.  Appendix 

A 

K Turbulent  kinetic  energy  per  unit  mass,  Eq.  (31) 

K Nondimensional  turbulent  kinetic  energy,  Eq.  (42) 

k Mass  fraction 

£ Length  of  inviscid  flow  boundary 

M Mach  number 
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m Mass  flow 

N Recompression  factor 

n Exponent  on  initial  boundary-layer  profile,  Eq.  (5) 

p Static  pressure 

Pt  Total  pressure 

R Nondimensional  radius,  r/rN 

Rg  Gas  constant 

Ri  Nondimensional  radius  to  intersection  of  inviscid  flow  boundaries 

Rt  Turbulent  Reynolds  number,  Eq.  (SO) 

r Radius  from  centerline  of  jet 

r„  Radius  to  intersection  of  inviscid  flow  boundaries 

T Static  temperature 

Tt  Total  temperature 

u Velocity 

X Nondimensional  distance  from  jet  exit  plane,  (x/rN ) 

Nondimensional  distance  to  intersection  of  inviscid  flow  boundaries 
Xi  Nondimensional  distance  to  beginning  of  recompression  zone 
X2  Nondimensional  distance  to  peak  pressure  in  recompression  zone 

x Distance 

x^  Distance  to  intersection  of  inviscid  flow  boundaries 
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xl  Distance  to  beginning  of  recompression  region 

X2  Distance  to  peak  pressure  in  recompression  zone 

x*  Distance  to  average  pressure  in  recompression  zone 

Y Mixing  zone  ordinate,  + to  high-speed  side 

y Initial  boundary-layer  ordinate 

a Angle  between  inviscid  flow  boundary  and  stagnation  streamline,  Eqs.  (82)  and 

(92) 

(3  Angle  between  inviscid  flow  boundary  at  (x^,  rj  relative  to  jet  centerline 

7 Ratio  of  specific  heats 

8 Boundary-layer  thickness 

8*  Boundary-layer  displacement  thickness,  Eq.  (28) 

e Eddy  viscosity 

r]  Nondimensional  mixing  zone  ordinate,  aY/x 

t?p  Mixing  zone  velocity  profile  shape  parameter,  Eq.  (4) 

i7b  Nondimensional  width  of  mixing  zone,  Eq.  (33) 

0 Boundary-layer  momentum  thickness,  Eq.  (29) 

p Density 

a Similarity  parameter  for  planar,  two-dimensional  turbulent  mixing 

cta  Similarity  parameter  for  axisymmetric  turbulent  mixing 

t Shear  stress 
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<(>  Nondimensional  velocity,  u/u^ 

0 Angle  between  stagnation  streamline  and  the  slip  line  at  xg,  Fig.  10 

00  Slip-line  angle  relative  to  jet  centeriine 
SUBSCRIPTS 

c.  Base  region 

D Dividing  streamline 

E Bleed  flow 

1 Initial  conditions 

L Low-speed  edge  of  mixing  zone 

m Conditions  in  mixing  zone  at  0 = 1/2 

N Jet  nozzle 

P Position 

S Stagnation  point 

U High-speed  edge  of  mixing  zone 

1 Jet  flow 

2 External  flow 

3 End  of  recompression  region 

- Inviscid  free  stream 


56 


